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I. INTRODUCTION 

The multichannel Kondo modelEI is the simplest im- 
purity model with non-Fermi liquid behavior. Originally 
introduced to describe "real metals" with magnetic im- 
purities, its applications go beyond the study of dilute 
magnetic alloys. For instance, it has been known for 
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some time that systems consisting of heavy atoms tun- 
neling between two neighboring sites and interacting with 
conduction electrons are a realization of the two-channel 
Kondo modelo. Another realization is the,quadrupolar 
Kondo effect in context of heavy fermionsEl. A detailed 
account of various aspects and applications of the multi- 
channel Kondo model is given in Ref. | 

For materials such as Pb\- x Ge x Te or K\- x Li x Cl al- 
loys, the tunneling may occur between an arbitrary num- 
ber of levels. Such systems could be modeled using a 
multichannel version of the Caablin-Schricffcr model, a 
SU(N) x SU(f) Kondo modclB. Here, N is the num- 
ber of spin degrees of freedom, and / is the number of 
channels, or flavor degrees of freedom. 

In this article we present an exact solution of the 
SU(N) x SU(f) Kondo model, and study the thermo- 
dynamic properties of the system. We obtain the leading 
exponents for the impurity contribution to the magnetic 
susceptibility and specific heat for arbitrary N and /. We 
also discuss the effects of channel anisotropy, which might 
drive the system from a fixed point with 7 = f /N > 1 
to a new fixed point where 7^ < 1. 

II. THE MODEL 

The Multichannel Coqblin-Schrieffer (MCCS) model 
describes electrons carrying two sets of internal degrees 
of freedom, to be denoted spin and flavor (or channel 
number) , interacting with an impurity carrying only spin. 
The impurity is localized at a point chosen to be the 
origin. The Hamiltonian reads, 

N f ,00 

W = -*EE/ ^l m (x)d x tp a . m (x) dx 

a m J -°° 
f N 2 -l N 

+ 2 E J ™ E E^U(°) ( T « D) ) Q >< m(0) 

m a. a.b 

dim(fj,) 

x E xl(W) a ,, b ,X»- ^ 

a' ,b' 

Both ipl m (x) and \ a are fermionic fields, the former cre- 
ates an electron at x with spin index a and flavor in- 
dex to, while the latter creates the impurity at x = 0. 
Imposing the condition ^2 a xi,Xa = 1, we have that 

Xa> [Ta^f Xb' represents the impurity spin operator 

V / a' ,b' 

in a representation of SU (N) specified by a particular 

choice of the matrices Ta &P \ where the index a runs 
from 1 to A 2 — 1, the number of generators of SU{N). We 
will restrict ourselves to the case in which the electrons 
are in the fundamental representation (denoted by (□)), 
and the impurity is in the totally symmetric representa- 
tion obtained from the direct product of /1 fundamental 
representations . 



The physical realizations discussed in the present pa- 
per correspond to the simplest case, [i = 1; the investi- 
gation of the /i > 1 cases gives us important insight into 
the general structure of the model and allows, in partic- 
ular, for a comparison with the results obtained for the 
multichannel Kondo model with impurity spin S > 1/2. 

In most of this paper we will study the isotropic model, 
J rn = J, with the symmetry U(l) charge x SU(N) s P in x 
SU(fY lavor . We will also assume that the different flavor 
levels are equally populated, N e = fN . 

In what follows shall solve the complete model and 
among other things, study its low-energy physics. As is 
well known, the low energy behavior of a system can of- 
ten be described in terms of effective hamiltonians, that 
are simpler than the starting hamiltonian; these are usu- 
ally referred to as fixed points. We shall determine their 
properties from the exact solution. We shall find that the 
model possesses a variety of fixed points (or low energy 
regimes), whose nature depends on the symmetry struc- 
ture in the flavor sector and on the spin representation 
(//), generalizing jthcu^miliar N — 2 case (the multichan- 
nel Kondo modefl'BEMJ. As previously, we shall identify 
the mechanism underlying the appearance of these fixed 
points as dynamical fusion by which electrons form spin 
complexes whose interaction with the impurity leads to 
a new behavior in the infrarecu. Each complex consists 
of / electrons fused into a local objects that transforms 
according to one-row Young Tableaux of length /. 

Within the Bethe-Ansatz approach a precise descrip- 
tion of the formation of these composites can be given. 
The linearized hamiltonian propagates separately the 
charge-spin-flavor degrees of freedom that make up the 
electron. Therefore the effect of flavor on the spin de- 
grees of freedom is recovered only in the physical space. 
To follow the dynamic coupling of spin and flavor we add 
some curvature which maintains the identity of the elec- 
tron while allowing its components to interact. It has the 
form H A = j K J2aJ2 m lZc 1 fa,^(x)dl'ipaMx)dx, where 
A is the curvature scale which is sent to infinity at the 
end of the calculation. Adding this term allows for the 
formation of bound states in the flavor singlet channel, 
which interact strongly with the impurity, and determine 
the low-energy dynamics even after the curvature is re- 
moved. A close analogy is a small magnetic field intro- 
duced to probe for magnetization, which may survive af- 
ter the field is removed. Imposing a cut-off D on the 
momentum variables guarantees the finiteness of the en- 
ergy. Other terms need to be added to the hamiltonian 
to maintain integrability, terms which we shall see below 
are irrelevant. 

Already for free fields the resulting theory is quite m- 
volved, and even the counting of states is not trivialQ. 
Nevertheless, the charge-spin-flavor separated basis is the 
natural one for the non-interacting problem, as we shall 
see later: it is the form to which the eigenstates tend 
when the interaction is turned off. We thus introduce 
the following elements: 
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A second derivative term with a curvature scale, A, 

1 f°° 
^A = ^EE ^l m (x)d 2 x ^ a>m (x)dx (2.2) 
Z1V a m J -°° 

which breaks charge-spin-flavor (CSF) separation 
of the linear spectrum. Once the electron compos- 
ites are formed, and the low-energy spectrum of the 
theory is identified, the scale is taken to infinity. 



JV -1 



Adding the term (2.2) also imposes restrictions 
on the form of the eigenstates which can be ex- 
pressed in terms of the following counterterms with- 
out which the model is not integrable for finite A: 

• An electron-electron interaction term, of the form 

/CO 

: I i III- i i li < ~~*~ > 



(2.3) 

When no impurity is present J can be chosen ar- 
bitrarily since the term has no effect on the linear 
spectrum. The linearized spectrum has a large de- 
generacy, and the inclusion of ( p^ ) and fl2.3| ) will 
provide a way to find the eigenstates. 

• A counterterm H CCl of the form 
^ = -lEE ^l m (x)V{x)^ a , m (x)dx (2.4) 



with 



V(x) 



■(S'(x +0 ) + S'(x-°)), 



(2.5) 



needs to be added to the Hamiltonian in order to 
preserve integrability at the origin; this term van- 
ishes once the curvature is removed, and plays no 
further role in the problem. 



A. First Quantized Hamiltonian 

A general Fock state of N e electrons and one impurity 
can be written in the following form: 

\f> = j: £ f \Rd Xj )Fi:^ b ({ Xl }) X i(o) 

{m j }{a j },b J -°° j 

In order for it to be an eigenstate the amplitudes F must 
satisfy the equation hF = EF, where the differential op- 
erator h, known as the first quantized form of the Hamil- 
tonian, takes the form 



j=l a 

N ° 1 

+ £ 2JS(xi - Xj )(Plj - V jt ) - J V ( x j), 

Kj 3=1 

with Pji(Vji) the spin(flavor) exchange operator, 

Pab,cd = S a dSbc, 
^mim2,m3m4 — &m\in^m>zm?,- 

The fundamental representation (□) is carried by the 
electron j and the (ji) representation by the impurity. 
When the latter is also in the fundamental representa- 
tion i.e. /i = 1 the hamiltonian can be rewritten as, 

h = J2 ~ id i + ( A_1 ) 9 j + 2J${xj)Pjo 

+ zj^i - x j)(Pij Tn) + £ j_ v{Xj) - 



Kj 



3=1 



B. S-matrices 

We will assume for now that both the electrons and 
the impurity are in the fundamental representation of 
SU(N). The eigenstate amplitudes are combinations of 
plane waves with pseudo-momenta kj, (j — l,...,N e ), 
and have coefficients that depend on the ordering of the 
electrons, and on the spin and the flavor indices. These 
coefficients are related through products of electron- 
impurity and electron-electron S-matrices that we will 
derive now. Consider first the wavefunction describing 
one electron (denote it by j) interacting with the impu- 
rity (denote it by 0), 



F 



™U i x j) = <? hiXj (KUA- X j) + B ZUA x j)) ■ ( 2 - 6 ) 



Applying h to it we have (we drop the indices in the 
amplitudes) 



k 2 

hF( Xj ) = (kj - ^)F( Xj ) 



(2.7) 



H(l - j-)(B -A) + JP j0 (B + A))5(x s ) 
-^(B - A)5'( Xj )e ik ^ + jV(x 3 )F(x 3 ). 

F is an eigenstate of h, with eigenvalue Ej — fc, ( l — 5^), 
if the terms in the second and third lines in (2.7) vanish. 
The last two terms cancel each other d ue t o the form of 
T5l) . The terms in the second line of (2.7) cancel if the 



amplitudes A and B are related by the electron-impurity 
S-matrix B = SjqA, where SjQ — (Sjo)^'"" is given by 
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i(l - kj/A) + JP j0 
kj/A) - JP j0 

i(l - kj/A) + J\ fi(l - kj/A) + JPj 



(2.8) 



i(l - kj/A) -Jj\ kj/A) + J 



Defining 



2J , v _ 1-x ( J 2 
c=- ~ , g(x) = - — II 



l-J 2 
we can write 



1 - J 2 V (1 - x) 2 J ' 



SjQ _ e -"^^ n ^rfm (a{kj/A) -icPjp 



g(kj/A) — ic ) ' 



(notice that arctan c = 2 arctan J) . Eventually we will 
send the cutoff to infinity. Therefore, expanding g(k/A) 
to first order in f/A 



l(k/A) ~ f - 



1 + J 2 \ fc 



l-J 2 A' 



we have 



Sj ~ e 



arctan / — 1 H~ icPj[ 



Xj — 1 + ic 



(2.9) 



where 



_ ' 1 + J 2 \ fc, 

- 1 w 2 " J a 



In the scaling limit, J and c have the same scaling be- 
havior. 

We now consider the case of two electrons. We gen- 
eralize the procedure followed in the case on one elec- 
tron: divide the configuration space into regions inside 
each of which there is no interaction and the wavefunc- 
tion is a superpositions of plane waves. They are six 
such regions in this case, corresponding to the ordering 
of three objects: two electrons and an impurity, and we 
label them by permutations Q e S3. For example, the 
element Q = (1,0,2) labels the region where electron 
1 is to the left of the impurity and electron 2 is to its 
right. We also introduce the notation 9(xq) to denote a 
function that takes the value 1 in the region Q and zero 
elsewhere. 

The two electron wave function is then of the form 
(Bethe-Ansatz), 

F?(x) = Ae^+ k ^ Yl d ( x Q) A a.,m, 

Q 

where m = (mi, TO2) and a = (0,1,02, do) and A is the 
antisymmetrizer. The amplitudes in the various regions 
are connected by S'-matrices, e.g. S 01 A 012 — A 102 , where 
S 01 , the electron- impurity S'-matrix has been already de- 
termined in the one electron problem. For this Ansatz to 
be consistent it must satisfy the Yang-Baxter Relations, 



gijgiOgjO _ gjOgiOgij 



(2.10) 



guaranteeing that the two paths from (1, 2, 0) to (0, 2, 1) 
yield the same answer. 

What is the electron-electron S'-matrix, S^l There is 
no direct electron-electron interaction term in the Hamil- 
tonian ( |2 . 1[ ) , and one may be tempted to adopt the naive 
choice S v — I for the scattering matrix of electrons i 
and j. Nevertheless, electron correlations are induced 
through the impurity. These show up immediately, since 
the naive choice does not satisfy the Yang-Baxter Rela- 
tions: and S i0 do not commute. This non commu- 
tativity captures some important aspects of the model: 
after electron i crossed the impurity the latter is left in a 
different state than before. Hence the state in which elec- 
tron j finds the impurity depends on whether it crosses 
the impurity before or after electron i. Herein lies the 
difference between a system of electrons interacting with 
a fixed potential (a one-body problem since all electron 
see the same potential) and a Kondo system, where the 
impurity correlates the motion of all electrons. 

Are we allowed to introduce a scattering matrix 
to satisfy the Yang-Baxter Relations? We proceed 
now to show that this is indeed the case, namely, the 
introduction of an electron-electron scattering matrix 
would not modify the original problem we set out to 
solve. Consider first the space of free electrons with 
a linearized Hamiltonian. The space is highly degen- 
erate: for example the energy E = k± + in the 
two electron space corresponds to a wave function F = 
j2 qe i(ki+<i)^+Hk2-q)x2 Aq for any choice of coefficients 
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basis of the form 



q . Equivalently, we can pick a 
F = e iklXl+ik2X2 (0(x\ - x 2 ) + S 12 8(x 2 - Xi)) A. The 
choice of S 12 is arbitrary in the two electron space, but if 
we wish to proceed to construct three (and more) electron 
wavefunctions then the scattering matrices satisfy must 
satisfy the YBE for electrons, S l3 S lk S jk = S^ k S ik S ij . 
When the Kondo interaction is turned on, the matrix 
is fixed by the interaction, which in turn picks the 
electron basis through the Yang-Baxter Relations (2.10). 

When the cut-off is present part of the degeneracy is 
removed already at the free electron level but the pro- 
cedure still goes through. Consider the model for two 
electrons away from the impurity. 



i 



1 



h = - idj -id l + - d '. + - dl 



2A 

2JS(x l -x j )(P lj -Vij) 



(2.11) 



This cut-off Hamiltonian is in the same universality class 
as the free linearized hamiltonian and possesses the same 
spectrum when the cut-off is sent to infinity; its particu- 
lar form was chosen so that the S-matrix it defines does 



indeed satisfy the (2.1C). Again, we divide configuration 
space into two regions: 

Skm+klXl ) (A^6{ Xl - x } ) + /C,"; 1 ^, - ««)), 
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and study the eigenvalue equation hF = EF. We have 



C. Eigenvalue equations 



hF 




(2-12) 

+[-i{B-A) -i(A-B) 
+i(k - *L)(fl - A))8(xj - Xl y^ +k '^ 
Vji){A + B)S( Xj - xi)e i to +kl * B * 

xi)e i( - kjXj+klX '\ 



1 ■((B-A) + (A-B))6'(x j 



2A 

The last line is identically zero; counterterms of the form 
( |2.4j ) are only necessary when the particles involved have 
different velocities. The rest of the terms proportional to 
S(xj —x{) cancel if the amplitudes in the different regions 
are related by the following electron-electron S-matrix, 



iaji + J(Pji - V 



jl) 



iotji - J{Pji-V. 



where ctji = (ki 
as 



kj) I A. Such a S-matrix can be written 



a 3 i 



2iJP, 



ji a ji 



2iJV 



ji 



Choosing 



Otjl 



J = 



2U 



J 



Otjl 



+ 2iJ 



1 + J 2 ' 

allows us to express the ^-matrix as, 



S 



ji 



Xj — Xi + icPji Xj — Xi — icP 



ji 



Xj — Xi + ic 



Xj - A; 



(2-13) 



(2.14) 



The S -matr ices ( ^!s| , 2.14 ) satisfy the Yang-Baxter condi- 
tions (2.10), and also 

gij gik Cjjk 

assuring that we were able to generate a cut-off version 
of the Hamiltonian while maintaining integrability. 

The cut-off scheme we introduced generates a flavor 
component in the electron-electron S-matrix. Clearly 
it captures the interaction among electrons induced by 
the impurity. More so, already for the free hamilto- 
nian Hq = -iJ2a Ei/-°l^L( :c )W«,mW a non 
trivial S-matrix must be introduced if we choose an 
SU(N) x SU(f) invariant basis (which is appropriate 
for a subsequent inclusion of an impurity interaction) 
rather than the simpler SlflfN). A careful counting 
of states can be carried outQ to show that all expected 
states then appear with the correct degeneracies. It is 
instructive that this would not be the case for the naive 
choice S/l, , = I. 

(flavor) 

The energy eigenvalues of a 7 V e -ele ctron state are a 
generalization of the first line of ( 2.12 ). They are of the 
form 



E 



N' , 



(2.15) 



In order to determine the spectrum, we impose pe- 
riodic boundary conditions, and solve the correspond-, 
ing eigenvalue problem. The procedure is standardly 
and we skip here the details. The result is contained in 
the Bethe Ansatz Equations (BAE) which we proceed to 
write down. Each of the degrees of freedom - charge, spin 
and flavor - is described by a set of variables whose num- 
ber depends on the symmetry of the particular state. The 
charge degrees of freedom are given by the set {kj,j — 
1, N e }. The spin degrees of freedom are parameterized 
by the sets {x^, 7=1, M r \ r = 1, N; M N = 0}. Fi- 
nally, the flavor degrees of freedom are represented by the 
sets {0^,7 = l,...,AP;r = 1,...,/;M/ = 0}. The set of 
integers M r ; r = 1, N — 1 specify the symmetry of the 
spin component of the wave function given by a SU(N) 
Young tableau with the length l r of the rth row given by 
V = M r - M r+1 , M N = 0, M° = N e + 1. Similarly, 
the quantum numbers {M r } specify the symmetry of the 
flavor component. 

The equations are 



M 

n 



x l i -{l-Xj) + U 



M l 

oj 1 



1 A7 



(1 - A,.) 



A; 



i 2 ;=i W 7 







— ic 


2,. 




1, 


<4 




+ ic 


<4 




— ic 


<4 


-ujj 


+ q 






- i- 




~ A/3 


+ ic 


X'y 


-x r 


— ic 


2,- 


..,7V- 


-1, 


X -y 


-Xf3 


+ ic 


X. ry 


-4 


— ic 


X.-v 




+ tf 




~x% 





n n 

t=r±l 0=1 



N e i 
3 = 1 7 



to' 



n n 

t=r±10=l 



x\ 



x\ 



e ikjL _ 

n 

0=1 
r = 

M 1 

n 

0=1 

M 2 

><n 

M r 

n 

0=1 

r = 
M 1 

n 

0=1 

M 2 

><n 

0=1 



The next step is to solve the equations for all possible 
states, identify the ground state and the low energy ex- 
citations. Subsequently, by summing over all excitation 
energies we shall obtain the partition function. 

The BAE contain the cutoff A which eventually is sent 
to infinity. We shall find that in this limit the equations 
reduce to a smaller set once the correct ground state has 
been identified. It is composed of string-solutions (see be- 
low) corresponding to electron composites which interact 



1 , ■ c N" 1 
Ai + *2 TT X, 



{l-\j) + i\ 



% i J=1 A n 



(1 - A,-) - t| 
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most efficiently with the impurity. To sharpen our intu- 
ition we begin by some strong coupling considerations. 



1. Casimin 



As mentioned, the mechanism underlying the physics 
of the multichannel Kondo model is the dynamic forma- 
tion of electron composites. We expect that those con- 
figurations are favored which allow minimization of the 
local interaction at the impurity site. Consider then the 
general problem of finding the ground state of the follow- 
ing hamiltonian: 



iv- 



j E T « )T « 



(2.16) 



where the set {T^, a = 1, N 2 — 1} is an arbitrary rep- 
resentation of SU(N), and {T*, a = 1, N 2 - 1} is the 
particular representation of the impurity, in our case it 
will typically be (/z). In this article we will consider im- 
purities with spin in a totally symmetric representation 
(for more general representations seetS) . Each set is nor- 
malized: Tr(T a T b ) = \5 a ^ b . 

The largest number of electrons allowed at the origin 
by the exclusion principle is N x /. This is obtained 
by placing N electrons in each of the channels. How- 
ever, such a state is a singlet, both in spin and in fla- 
vor, and gives a zero contribution to (2.16). Therefore, 
the number of electrons that form the composite, M, is 
such that M < (N - 1) X /. We will show here that if 
the impurity is in a totally symmetri c rep resentation, the 
electron composite that minimizes ( 2.16 ) is made out of 
M = (N - 1) x / electrons. 

We will characterize the different representations of 
SU(N) by their Young tableaux. The fundamental rep- 
resentation is denoted by a box, □, and the singlet by a 
point •. The totally (anti)symmetric representation re- 
sulting from the direct product of [i representations is 
denoted by a single (column)row made out of boxes, 
where in the antisymmetric case we assume /i < N. 



Antisymmetric 
t 

I □ 



Symmetric 



mi 
m 2 



When N = 2, the interaction (2.16) can be written in 
terms of conserved quantities 

JS Me ■ S l = 

The operator S(S + 1) is a particular case of the Casimir 
operator, C(T), of SU(N), which commutes with all the 
generators of the group. For arbitrary N, 

jJ2 T ™ T a = i ( c ( rT ) ~ c( rM ) - c (n) ■ (2.18) 



Given a representation T of SU(N) with M-.boxes dis- 
tributed according to the set {mj}, we haveEj 



C(T) 



M(N 2 - M) 
2N 



1 k 



+ l-2j). (2.19) 



i=i 



We can use Young tableaux as an easy way of decom- 
posing the direct product of representations into a-direct 
sum. The procedure is standard (see for instancecJ). 

The electrons will be evenly distributed among the 
channels, forming a flavor singlet. Hence, the spin of the 
electron composites is described by a rectangular tableau 
with / columns and k < N — 1 rows. Multiplying the 
electron tableau by the impurity tableau we have (in the 
graphic representation we drop from the tableaux the 
singlet part consisting of columns of length N), 



f A 4 
rrm 



fc + 1 



k + l 



n > f 



An arbitrary representation resulting from a product of 
M fundamentals is associated with a Young tableau made 
up of M boxes, distributed in k < N rows. Let rrij be 
the number of boxes in the j-th row. Then, rrij > m J+ i, 

y^,j—i mj — M ■ The corresponding Young tableau will 
be of the form 



fc + 1 



A» < / 



for k < N — 1. If k = N - 1, we have 
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N - i ^flm, ® mm 



cm + 



N - 1 



A» > / 
/* = / 

A» < / 



Notice that we have drawn only the terms in the de- 
compositions that give the lowest energy. The energy for 
all such configurations is given by 

jJ2T a M ^ = -k^^Il(N + m a X ^J)). (2.20) 

a 

Therefore, the energy is minimized when composites of 
(N — 1) x / electrons are formed. 

There are three different situations depending on the 
value of n/f, as in the multichannel Kondo problemEl. 
When fj, > /, there is underscreening: the electrons can- 
not screen the impurity completely and the spin con- 
figuration is characterized by a Young tableau with one 
row and [i — f columns. As in the N — 2 case, we will 
see later that such object behaves as a free spin in the 
Kondo problem. The second case, fi — /, corresponds 
to complete screening: the electrons and the impurity 
form a singlet. This is a stable fixed point of the full 
Hamiltonian with Fermi liquid behavior. Finally, fi > f 
corresponds to overscreening: there are more electrons 
than necessary to screen the impurity. The resulting ob- 
ject corresponds to a tableau with N — 1 columns and 
/ — p rows. This configuration is unstable to the kinetic 
term, and the fixed point in this case is characterized by 
Non-Fermi liquid behavior, as we will see later. 



2. Fusion 

We turn now to the dynamics of the full model cap- 
tured by the BAE. We shall argue that the ground state 
and low lying excitations lie in a sector of the theory 
given by solutions of a particular form - / — strings. So- 
lutions of this type are SU(f) flavor singlets - allowing 
them to have maximally large SU (N) spin. We shall find 
that this class of excitations is characterized by a scale 
To = De~7^ . When strings are broken to form flavored 
excitations we expect them to be characterized by other 
scales which will tend to infinity as the cut-off is removed 
and thus not contribute to the impurity dynamical 

The formation of composites in flavor corresponds to 
solutions of the BAE where the charge parameters, 
are complex numbers centered around {w^,}, according to 
the string hypothcsisotl. Likewise, rank r flavor parame- 
ters are themselves centered around rank r+1 solutionsB. 
The form of the charge parameters is, 



q= 1,2,...,/, ps real. 



while the flavor parameters, 

7 = 1,2, M r } = { PA / A + U[(f - r + l)/2 - q]; 

q=l,2,...,.f-r,A=l,...,N} 

where r = 0,1,...,/ — 1. These configurations satisfy 
the BAE in a trivial manner and induce fusion in the 
BAE equations as well as in the form of the wavefunc- 
tions. A string built on momentum p as its real part 
induces in the wave function a composite of the form 
exp{— ~AJ ^ j \xj —xi\ + ip(x\ + ... + Xf)} x [...], which 
becomes local as A — > oo. 

Inserting the string configurations into the full BAE 
we obtain the effective equations governing the impurity 
spin dynamics. After removing the cutoff they become, 



JfpsL 



M i 

X 



TT 7 



1 + ifi 



7=1 X ~t 



M" r 
13=1 X ~< 



x' 



M* 



ic 



n 

t=r±l/3=l A T 



t ■ c ' 



X + ic _ Xj + ■ 



r = 2,...,N-l, 



M 1 1 

TT — 



TT 

11 yi-l 



i+in 



M yl 

><n" 
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xl + H 



13=1 X ~< 



1 ■ c ' 
Xp-H 



and the energy is given by 



E 



N e /f 



5=1 



We proceed now to discuss the solutions of the fused 
equations. The solutions for the rank r spin variables 
{x r } again fall into strings of arbitrary length n, 



X, 



X r ~i n + + 1 - 2j); j = 1, n, n = l, oo. 



and a state is characterized by the quantum numbers 
M r ' m specifying the number of length-m strings of rank 

rATZ=i mMr ' m = M [)- 

The equations coupling the real part of the strings, 
after summing over the complex variables, can be conve- 
niently written down in a logarithmic form. Let us first 
introduce the following definitions: 



6 n (x) 

4>n,m( x ) = 



-2 arctan | — x 
nc 



min(n,m) 

E • 

min(n— l.m— 1) 



i+n+k-2j{x), 



3=1 



^m+n— 2j (*^) • 



7 



Then, after some manipulations the Bethe-Ansatz equa- 
tions take the form, 



oo M r ' m 

E E (< m (xr-x r n+€, m (xr-x r e m ) 

m=l /3=1 

oo M r_1 ' ! 

2 ^"+E E <t>Uxr-x r f u ) 

1 = 1 (3=1 

oo M r+1,! 

E E <i(xr-x r f3 +1 ' l i 

1=1 (3=1 
oo M 1,m 

E E (<m(34 ,n -xj ,m )+< ro 04' n -xj ,fn ) 

m=l /3=1 

oo M 2 ' 1 
i=i /3=i 



The expression for the energy of the spin and charge 
sector is given by 

EZtt 
-yms 



5=1 



D 



oo M 1,n 



7 E E «/(x^ l -l)-^min(n,/)) 
n=l 0=1 

where D = ±j- is the electron density. It will turn out 
also to play the role of the cut-off. In the presence of a 
magnetic field, H, there is a contribution to the energy 
of the form, 



w_1 N-^ 
-2H E (M k - M k+1 )(— k) = 



k=0 



-2H 



N - 1 N 

(^-(^v e + i)-E Mr 



We now take the thermodynamic limit, N e — > oo, 
L — > oo, holding D finite. In the limit we may replace 
sums with integrals after introducing densities of solu- 
tions, cr^(x), and densities of holes in the distribution of 
solutions, cr£> (x)- The energy is now written as 

E = E c - H(N - l)(N e + 1) 

oo N-l 



^ ^ /*oo 

EE/ d XVn(x)9rAx), 



and -Ec denotes the contribution of the charge sector to 
the energy, 

N "/f 2n 
E c = E ~L m& - 

(5=1 

In the thermodynamic limit the BAE are replaced by 
integral equations for the densities, {a r n , a% h }. Standard 
manipulations!!^ lead to, 

oo N 

--EE/ d xX'f m (x - x'Xix')- 

m =ls=0 J -°° 

where we introduced the following operators: 

min(n,m) 

K n,rn(x) = E k ^+n+ a -2j (x), 

min(n — l,m— 1) 
^°,m(x) = E k m+n-2j(x), 

1 



Mx) = -^ a (x), k (x)=6( x ). 



and 



B n ,m{x) = K^m(x) 

A r n s m = A n<m 6 r < s - B n ^ s+1 + S^- 1 ), r > 1, 



and by convention, 



N e 



where we introduced the energy function, 
3rM) = J «/(X - 1) - 7T min(n, /)) S rA + 2Hn, and itg contribution to the spin free energy 



<r°(x) = <5( X )<5 U + 6( X + 1)*,,,, (2.21) 

^( X ) = 0. (2.22) 

We shall not analyze here the ground state and the 
individual excitations. Instead, we shall proceed to derive 
the thermodynamic properties of the model. 



III. THERMODYNAMICS 

A. Thermodynamic Bethe Ansatz equations 

We now calculate the impurity contribution to the free 
energiL. using the well known formalism of Yang and 
Yanguu. We seek to find the configuration {c^(x) + 
°ti \X)} which would extremize the free energy. The 
entropy of such a configuration is, 

S = Y,j d x{ KM + < h (x)) In «(X) + < Jl (x)) 

n.r 

-<(x)ln<(x)-< Jl (x)lna^( X )}, 
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F = E — TS = 

-H(N-l)(N a + l) + J2 [ d X {<(x)9rAx) 

n,r 

-T + a^ix)) In «(x) + <^(x)) 

-< (x) In < (x) - at (X) In (*))}. 

The free energy is varied with respect to the densi- 
ties, subject to constraints imposed by the Bethe-Ansatz 
equations, 



S*1(x) = 6<t»(x) = 0. 

We obtain the following infinite set of integral equations 
for the equilibrium densities, 

ln(l + r ? ;(x))= (3.1) 

9r,n(X) 



T 



where 



/oo 
dx'A^ix-x^Hi + ivKx'))- 1 ) 
... . -oo 



rfnix) 



_ a ri (X) , , .V . I. _ 11 , I _ 



<{x) 



, to-* = = o. 



We transform this set of equations with the help of the 
following identities (we will now drop the functional de- 
pendence) 

Ar,s _ n( A r - S i A"^ s \ 

n,m \ n—X,m n+l,m/ 

= S ntm S r ' s - GS n , m (5 r ' s+1 + S r ' s ~ 1 ) 

A r C - GA r 2 > s rn = 5 hm 5 r > s - G5i, m (5 r > s+1 + J^" 1 ), 

with the integral operator G defined as 

[i] _ i r , , fix 1 ) 



Gfix) 



i r°° 
- f(x) = - / d x '- 

\ 2c J-oo < 



a refund 

-Hi + ( V ; 1 )- 1 ) = f-s n j 

N-l 

+ J2G r N Q (Ml + ri q n + l)+Ml+V q n-l) 
9=1 

-G-Hnil+vD), (3.5) 

where ln(l + t]q) = 0, and the Fourier transform of the 
kernel of the integral operator G 7 ^? is given by 



s< h (x) = - E r d ^' A ^mbc - x')s<(x'l g n(p) = 

™ „ J — oo 



_ sinh (min(r, q)^f) sinh ((iV - max(r, q))^) 



[0] + [2] J w ~ 2c J_ x ^ cosh (f ( X - X')) ' 
and [n]f( X ) = JZo d x'Mx - x')/(x')- We find, 
h "/i = -y^arctane^(x- 1 )^ 1 ,5„j + Gln(l + ^) 

-G(m(l + far 1 ) -1 ) + ln(l + i^ 1 )- 1 )), (3-2) 

+G(ln(l + < l _ 1 )+ln(l + < 1+1 )) 
-G(ln(l + [rfn 1 )- 1 ) + ln(l + (r?^ 1 )" 1 )), (3.3) 
with boundary conditions, 

Urn ([n + 1] ln(l + rQ - [n] ln(l + rf n+1 )) = -2^. (3.4) 



which follow directly from eqns ( |3.l[ ). Another form of 
these equation can be obtained after inverting, see Ref. [ 



sinh (f) sinh (JVf ) 
The driving term in these equations, 

A^; id = G^G- 1 (^jarctane^*- 1 ) 



is the energy of the fundamental excitation. It can be 
calculated explicitly, 

a refund D ( N -r 

( fir N — r\ . ( 7r . 
- 2 arctan I tan I — — ^— I tanh(^— (x - 



f-2 f-1 



a 



-o — — ^- 



FIG. 1. Diagrammatic representation of the integral equa- 
tions. The circles correspond to the functions r[ n . The filled 
circle indicates that the equation for the corresponding rf n 
has a driving term. The circles with stripes indicate that the 
corresponding rf^ have driving terms in the other set of TBA 
equations. The solid line indicates a link between two rf n 
thorough the convolution Gln(l + 77). The dashed line indi- 
cates a link through Gln(l + 1 / 77) . Finally, the box encircles 
the r/1, which are the functions used to evaluate the impurity 
contribution to the free energy. 



A pictorial description of (3.5) is shown in Fig. [I]. 
The circles correspond to the functions rf^ix), and are 
arranged according to their indices. The lines join func- 
tions that appear in the same equation. The full and 
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dotted lines indicate that t he functions f]^ 1 and T]n±i 
app ear diff erently in eqns (3.2, 3.5). The driving terms 
in (3.2 3~3| ) correspond to filled dark circles in Fig. [|. 
The diagonally lined circles cor resp ond to functions as- 
sociated to a driving term in ( |3.5| ). Clearly there are 
two regions: the one, n < f, contains a finite number 
of functions rf nl while the other, n > f, is unbounded. 
The regions are separated by the column with the driv- 
ing terms, n — /. When studying the low-T properties 
of the system we will only need to consider one region at 
a time. 

We will now write the free energy in terms of the set 
{t/£}. Using the integral equations for the densities we 
can write 

F = F + J2fd X (9r,n<r r n - T< ln(l + fQ 

n,r 

+ T E / tt^lnfl + ti)- 1 ) 

where F =E C - H(N - l)(N e + 1) is the ground state 
energy. After a few further manipulations the free energy 
can be written as, 

/oo 
dxMl + iVn)- 1 ) >< 
-°° 

N e min (™</) 



We are only interested in the impurity contribution to 
the free energy, F l , which contains all the effects of the 
interaction. It is 

OO />oo 

F' = -TJ2 / dxSixMHl + ivi)- 1 ) 

n=1 J-co 
N-l r oo 

= E / dxG^(x)gl( X ) 

q=l J -°° 
N-l r oc 

q=l J-oa 

The first term corresponds to the impurity contribution 
to the ground state. At finite temperatures we are only 
interested in the second term which after further manip- 
ulations becomes, 



N-l 



F* = -T £ d X 



1 sin7r 



N-q 
N 



^ J_ 00 Nc cosh ^ + cos 7T 



N-q 
N 



ln(l+^ 



When the impurity is in the fundamental representation, 
only the j]\ functions contribute to F\ In Fig. [I] this 
feature corresponds to a box drawn around the first col- 
umn. 



B. Scaling limit 

We will now take the scaling limit, D — ► oo, c — » 0, To 
constant, where 

T = De~^. 

This is the correct limit as discussed in Ref . [ || . We also 
introduce the new variable 

t 2tt T 
^ Nc T 



so 



F' 



2tt 



(smvrV)ln(l + ^(0) 
uH ,-oo cosh(^-ln^) + costt-^' 



(3.6) 



The only modification in the thermodynamic equations 
is in AE f ™ d . Thus 



HI + (.Vn) ) 



f 



sin 



(fK 



(3.7) 



JV-1 



+ J2 c^'CMi + €+i) + + €-i) 

-G-Hnil + vD), 
ri r =0, ln(l + (^r 1 ) = 0, ln(l + (^r 1 ) = 0, (3. 



with boundary conditions 



H 



hm ([n + 1] ln(l + rf ) - [n] ln(l + = -2-. (3.9) 

n— >oc J 



C. Asymptotic solutions-Low temperature properties 

Here we will study several asymptotic limits of the 
thermodynamic integral equations. Some technical 
points will be considered in detail. It is not e asy to study 
analytically the integral equations (3.7-3.9) due to the 



complexity of the operator Gj?_. Instead, we will study 



the equivalent set equations (3.2-3.4). We will discuss the 
appropriate procedure to obtain the asymptotic solutions 
of the equations order by order. 

The zeroth order approximation yields a description of 
the fixed point itself, the corrections (first order) describe 
its neighborhood. 



1. Zeroth order - the fixed point 

The functions if n tend either to or to constant val- 
ues as the magnitudes of their arguments tend to in- 
finity. The only information needed about the driving 
term is that it tends to as £ — > — oo, and to — oo as 
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£ — > oo. T herefore , the parameter x does not appear ex- 
plicitly in ( 3^- 3^ ) , and the kernel of G can be replaced 
by (l/2)«f(£- £'). 

Thus, the zeroth order problem consists in evaluating 
the following set of constants 



n, 



r,± 



t£K - ±oo). 



When £ — > — oo, all the driving terms vanish, and the 
algebraic equations for the set {rfc~} are 



21nT#- = ln(l + rfej + ln(l + ^Ci) 

+ ln(l + (^+ 1 -)- 1 ) + ln(l + « r 1 --)- 1 ), 

7 /o " = °> 

ln(l±«-)- 1 )=ln(l + (^-)- 1 ) = 0, 



) 

T 



lim ([n + 1] ln(l + <•") - [n] ln(l + ^)) 



Since the kernel of G^ (x) satisfies 



/~*r.q s~iN — r.N — q 



we must have 7y„ r — rf n . The solution is easily obtained 
as it is independent of the flavor symmetry, 



sinh((n + r)xo) sinh((n + N — r)xo) 
sinh(riro) sinh((iV — r)xo) 

n = l,2,..., r=l,...,N 



- 1, 



H 



1, x = - (3.10) 



As for 77™'+, we should consider separately the cases 
n > f and n < f (?7y + = 0). In the first case, we 
have equations similar to those for {?7^' _ }, except that 
all the indices n are shifted by /, in analogy with the 
multichannel Kondo model. Hence 

, sinh((n — / + r)x ) sinh((n - / + N - r)x ) 

Tj = — 1 , 

smh(rxo) smh((N — r)xo) 
n = /,/ + l,..., r=l,...,JV-l. (3.11) 

Finally, for n < f there are a finite number of in- 
volved, since rff + — 0, and Gln(l + iff) — 0. As in the 
multichannel Kondoi, the sink functions are replaced by 
sin functions, and the coefficients are independent of the 
magnetic field. Thus 



sHxjfN)( n + r )) sin ((7Tlv) (" + N ~ r )) 



1, 



sin( I7 ^ J r)sin( I7 ^ J (7V-r)) 
n = l,...,/-l r = l,...,JV-l. (3.12) 

These results coincide with those obtained in Ref . [ [l5| for 
a model of interacting fermions with the same symmetry. 
This is not surprising since these results depend only on 
the symmetry of the problem. No tice that t he multichan- 
nel Kondo results correspond to (3.1C-3.12) with N = 2. 
Some features of ( 3. 10| - 3 . IS ) can be appreciated in Fig. 
[l]. When £ -> -oo the driving term does not contribute 



to the equations and the situation is the same as in the 
Coqblin-Schrieffer model. For £ — > oo, and n > /, we can 
disregard the n < f sector, and the leftover diagram is 
effectively the same as for the Coqblin-Schrieffer model 
with the substitution n — > n — f . Finally, for n < f we 
have a finite number of nl; + involved. Hence the replace- 
ment of the sinh by sin. 



2. Residual entropy - the fixed point 

Here we will calculate the residual entropy in the over- 
screened case, / > 1. As T — » 0, the dominant term in 
the free energy will be linear, and it will depend only on 
the values of rj[ ,+ . 



rri N-l 

2tt ^ 



sin 7r 



q=l 



N -q 
N 

ln(l + ^ + ) 



cosh(£ - In 2i) + cos (tt^) 

jv-i 



-T 



o=l v 7 



ln(l+^+), 



(3.13) 
(3.14) 

(3.15) 



Substituting the values of n\ ,+ , and taking advantage of 
the symmetry that n 1 ^ q,+ = nf' + , we find 



F l = -Tin- 



ttJV 
f+N 



sm 



f+N 



Hence, the residual entropy is 



q i dF* 



= In- 



ttN 
f+N 



= In- 



ZLf 
f+N 



T=0 



f+N 



sm 



(3.16) 



f+N 



Once again, we recover the multichannel results if we set 
N = 2. It is quite clear that it is not the logarithm of an 
integer number! 

The expression for the entropy can be written as the 
sum of two terms: one that depends only on N + f and 
a second one that depends only on | log7|, (7 — f/N) 



S l = In sin- 



In sin ■ 



l + el lo s-H N + f 

In the limits / 3> N and /<JVwe have 



(3.17) 



P 



In/ 



f<ZN. 



(3.18) 



6 N 



Furthermore, it is clear from (3.17) that two systems 
characterized by 71 and 72 such that 71 = I/72 have 
the same residual entropy. When 7=1, the first term is 
zero. 
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Fig. H corresponds to (3.16) for different values of N 
and /. It is quite apparent that the value of S l increases 
with N + f. It is also clear that the figure is symmetric 
with respect to the N — f axis, which means that S' 1 
is the same for 7 and for 1/7. Finally, if we fix N + f, 
the largest value of the residual entropy corresponds to 
N = f. 




FIG. 2. Overscreening residual entropy, S 1 , for an impurity 
in the fundamental representation of SU(N), and for different 
values of N and /. 

In terms of the diagrammatic construction, Fig. [I], 
S l measures the size of the overscreened region and how 
asymmetric the region is. For fixed N + f, the largest 
residual entropy corresponds to 7 = 1, in the same way as 
the square is the rectangle with the largest area for a fixed 
perimeter. This can be seen in the following diagrams, 
were we have omitted the lines and drawn only the circles 
corresponding to n < f. 

N = 7,f = 3 

o o 

o o 

o o 

o o 

o o 

o o 

(3.19) 



N = 


/ = 


5 








O 


O 


O 


O 






O 


O 


O 


O 






O 


O 


O 


O 






O 


O 


O 


O 






N = 


3,/ 


= 7 








O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 


O 



(3.20) 



(3.21) 

The largest entropy corresponds to the configuration 
with the largest number of circles, for a fixed value of 
N + f. That is, N = f. Notice that the first and the 
third cases have the same number of circles and, indeed, 
the same value of of the impurity entropy S l . 



3. First order - the neighborhood of the fixed point 

Now we turn to the calculation of the thermodynamic 
properties of the SU(N) x SU(f) mod el we ll below the 
Kondo scale T . As is obvious from Eq. ( 3.14 ) for T«To 
the nontrivial temperature dependence of the impurity 
free energy is determined by the asymptotic behavior of 
the functions r/K^) in the region £ ^> 1. Therefore we 
want to find the dominant dependence of rf n on £ for 
n < f and £ large and positive. To this purpose we will 
only need the equations with n < /, which do not have 
a driving term, and the asymptotic value of iff. 

Consider the action of the operator G 



XJ 2c 
1 

~ 2c 



N 

47T 



1 



fix') 



cosh(f(x-x')) 
cosh(f^( X -x')) 



cosh(f (£-£')) 



where we have dropped the tilde in the last expression. 
This establishes the correspondence between functions 
and variables of the two sys tems of equations. 

Fro m the set of equations (|0-3.4), and the asymptotic 
values ( |3.7| - 3.9 ) we learn that for large and positive values 
of £ we have (if xq = H/T is very small) 
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t£(0 oc (a r 4 (3 r xl)e-7 si <^> 



(3.22) 



We now evaluate the dominant contribution to ln(l 
ft) 



/OO 



, Mi + ??p 

cosh f 



^(a p +/3 r xg) / d£'- 



cosh f e ' 



where w is some lower cutoff of the integral of the order 
of unity and 

For large £, the only relevant contribution to the integral 
occurs around ~ 1/A. Therefore, we approximate the 
previous integral by 



(a r +f3 r x 2 ) 



l + A N 



{a r + f3 r x 2 )e 



N c 

Notice that this is correct up to terms of the form ^e~~* 
which cannot be accounted for using this crude approxi- 
mation. 

The previous calculation indicates that ff n< j will have 

a contribution of order e - ^, since G^e^ a( - oc e~ Q ^, 
as we will see below. Therefore, we have to determine 
whether there are contributions more singular still. In 
other words, we have to find out if there are solutions of 
the integral equations for large £ of the form 



^</(0~<</ + <(a r + /3 r ^)e- T «, with r< 



N 

'f+ c 'n( a '+P' x o) e "*) With T < —. 

Introducing the eigenvalue A = 2 cos ^ , so that 
Ge~ T « 



A 

we proceed to convert the TBA equations into algebraic 
recursion relations. Noting that, 



ln(l + n r n ) ~ ln(l + ft+) 4 (b r n 4 a r n xl)e-^, 
\nr) r n ~ h£ + 4- 

where, 



•,+ , ( b n + a n X l) e -rj 



(b r n + a r n x 2 ) 



r,+ 



c r n (a r + P r x 2 ) 



l + i} 



then substituting in the integral equations (3.2,3.3) for 
n < f and using the zeroth-order results we obtain the 
following set of algebraic equations for the coefficients of 



e T ^ (we only write the equations for b r n since they are 
identical to those for a^), 



a&; = + b r n _ x ) + -j^b r n +1 + -^b 



with 



b r = b r f = 0. 



n s 



(3.23) 



More explicitly, upon inserting zero order values the 
equations become, 



Xbl 



sin((n + N)a) sin(na) 



sin((n 4 


r)a) sin((n - 




-r)a) 


sin((r 4 


l)a) sin((7V 


— r 


-l)a) 


sin((n 4 


r)a) sin((n - 


~N 


-r)a) 


sin((r — 


l)a) sin((7V 


— r 


4l)a) 


sin((n 4 


r)a) sin((n - 




-r)a) 



b r+l 



(3.24) 



where 



f + N 



We solve ( 3.24 ) by ins pecti on. Since b r n has to satisfy the 
boundary conditions, ( |3.23 ), we have that 



b r n = sin((n 4 N)a) sin(na)rf^ 

is the maximal solution when d r n = d = constant, and 
the eigenvalue is, 

2tt 



ITT 

A = 2cos- = 2cos / + iV . 



Hence, finally 



2N 



(3.25) 



(3.26) 



N + f 

We shall see next section that t is the main critical ex- 
ponent in the model. 

4- Specific heat and finite temperature susceptibility 

The expression for the impurity contribution to the 
free energy, F l 7 has always at low temperatures a term 
which is proportional to T 2 . 

This contribution comes from the term proportional 
to e~~^ present in n r n< j, as we discussed in the previous 
subsection. 

Here, we will study contributions of the form 



rn poo n sin 
AF Z oc -— / d^ 



J v ^ cosh (£ - In 2i) 4 cos ir^ 



rp / rp \ T pOO 

27 \Y 



r,oo ^sin(Tr^) (a* + /3»z§)< 



(3.27) 
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which might become dominant depending on the value of r. We will consider the three cases f > N , f = N , and 
f < N separately. 

Case f > N: In this case, t < 1, and 



lim 

5^-00 cosh £ 



0. 



That means that we can make the same the approximation for the free energy that we made when we evaluated 
G(: ln(l + rf^). Therefore, we have 



- 'TV 

n, 



(3.28) 



T l TS l - T(A + B 

with A and 5 being constants o the order of unity and we obtain 

N-f 

(3.29) 

Needless to say, wi^p^Tf^— 2 we recover the multichannel results. As a matter of fact, the exponents depend only on 
the ratio 7 = f/N 



2N 
J 1 \ TT+7 



x x tAt 



Case f = N : Since r = 1 in this case, we cannot extend the integral in eq. (3.27) to £ — > —00 and we have to restrict 
it to the interval [S, 00), where S is a finite number of the order of one. Making use of J jj^zt = z — \ ln(l + e 2z ) we 
have for very low temperatures 



TIT 
AF l oc — 



N 



(sin»£^)(a«+/)'i»)e- 

N-q 



2n \T JJ S _i n 2L f-J cosh (£) + COSTr-^ 2 



-TS 1 



T 2 InT 



(A + B 



Hence 



T T IT 

C oc - — In — , y 1 oc In — . 

T T ' T To 



(3.30) 



Case f < N: In this region 1 < t < 2. 
Consider the integral 



1 



5-i og lo l + e 2x t 1 
2 



(/,!/ 



(3.31) 



It is possible to find a primitive for —^j integer; we have, 



f & £2=1 cos 7f(2 "if fe ' 1) Ml - 2y cos ^ + v 2 ) 



dy 



y 2 (l +y— 1) 



1 

-- + < 



7r(2fc-l)(2w-l) 
2n 



arctan 



y — cos 



2n 

Tf2fc-1) 



r-1 



2n 



ln(l+y) 
2n+l 1 2n+ 



TT £fc=i - ^K^l) ln(1 _ 2y cos + ft 



7r(2fc-l)(2w) 



2n+l X)fc=l Sm ' 2n+l 

As T — > 0, the leading terms are of the form 



2n+l 

arctan 



ir(2fc-l) 
^~ COS 2n + l 

. 7T(2fc-l) 



= 2n+ 1 



T-l 



ctn. 
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T 2 / T 

— cin. T — 

T \n 



Since r > 1, the dominant term in AF Z is of order T . 
As for the specific heat and susceptibility, we have 



T ( T 



2N 
N + f 



1 / T 



N-f 



(3.32) 



These results are valid for any 1 < t < 2, as can be 



verified by numerical integration of (3.31), or by the nu- 
merical solution of the thermodynamic equations. 

To summarize, there are three different kinds of behav- 
ior in the overscreened sector, depending on the value of 
the ratio 7 = f/N. i) When 7 > 1, both C l /T and x i 
have power-law divergences as T — > 0. The behavior is 
similar to that of the multichannel Kondo model with 
/ > 2. Indeed, the exponents are the same, since they 
depend on 7 only, ii) For 7 = 1 there are logarithmic di- 
vergences as in the two channel Kondo model, iii) When 
7 < 1, the values of C' l /T and \ % at T — are finite. 
Actually, it can be deduced from the numerical analysis 
that these constants are the same as in the corresponding 
completely screened cases (// = /, as we will see later). 
However, the fixed point has Non-Fermi liquid behavior, 
as can be seen from the value of the residual entropy and 
from the subleading power-law terms. 

One can relate the different kinds of behavior to the 
shape of t he n < f s ector in Fig. [l] as can be seen in the 
diagrams ( 3Tjj- 3.21 ) . The square diagram corresponds to 
7 = 1, whereas the horizontal(vertical) one corresponds 



D. Channel Anisotropy 

In this section we consider briefly the case when some 
of the couplings J m are different. From the study of the 
analogous problem in the multichannel Kondo modeO, 
we conclude that up to / different energy scales will ap- 
pear in the problem depending on the pattern of sym- 
metry breaking. The novelty here is that there might 
be a situation where 7 > 1 for an intermediate regime 
of temperatures, whereas for very low temperatures the 
behavior is characterized by an effective 7 smaller than 
1. 

Consider a system where the flavor symmetry is such 
that p energy scales T\ < T2 < ... < T p , are generated. 
Each scale Tj is related to a driving term at the level 



in the TBA equations (3.7). We will assume for 



simplicity, that mi < 7712 < ... < rri p = /. If the largest 
flavor symmetry possible is SU(f) there will always be a 
driving term at the level n = f . 

Then, when the temperature is belo w any Tj, the ther- 



modynamic properties are given by ( pug ) and ( j3~29l ), 
where 7 = f/N is replaced by j e ff = mx/N. As the tem- 
perature is increased, the behavior of the system when 



Tj-i < T < Tj corresponds to 7 e // = rtij/N. Indeed, the 
value to the impurity contribution to the entropy will be 
close to S l (rrij , N) . 

Flavor anisotropy is a relevant perturbation of the 
isotropic hamiltonian. In general, the system will flow 
away from the fixed point characterized by /, and N to 
a, new fixed point characterized by mi < /, and N. From 
( 0.16 ) we see that <S* is reduced in such flow (S l (f, N) is 
monotonous in both / and N). 

It is worth noticing that, once mi < N, \ % an d C l /T 
become constant as T — » 0. 






FIG. 3. Same as Fig. 1, but now the interaction with the 
impurity breaks the flavor symmetry from SU(f) down to 



Y[ P =1 SU(irij — rrij-i), with mF- = /, 



m 



0. 



The system of TBA equations are represented dia- 
grammatically in Fig. ^|. Depending on the pattern of 
symmetry breaking, driving terms appear at different val- 
ues of n. There are always driving terms for n = f. The 
properties a temperatures between two different scales 
are related to the corresponding overscreened part of the 
diagram. Notice that channel anisotropy may have the 
effect of changing the s hape in the overscreene d area from 
somet hing similar to ( 3.21 ) to diagrams like ( 3.20 ), and 
( |3.2lD , but not the other way around. 



E. Impurity in a higher dimensional representation 

Finally, we study a generalization of the model in 
which the impurity behavior is that of an object in a 
rank /i representation of SU(N). In the SU(2) case it 
corresponds to an impurity with spin S. Following Ref. [ 
plU and the same formulation that we followed for the 
fundamental representation, we find the following set of 
effective Bethe Ansatz equations 



JfpsL 



M 1 1 



l + ifi 



M r r 



n n^f- 

r = 2,...,A-l, 
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FIG. 4. Same as Fig. 1, but now the impurity contribution 
to the free energy involves the set ryJJ, where n is the rank. 

The impurity contribution to the free energy is 

/oc 

As in the /i = 1, a series of transformations allow us to 
rewrite the free energy as 



F l = 



JV-1 

E 

9=1 



9=1 ^ X 



d x G^( X )ln(l + <) 



where the first term corresponds to the impurity contri- 
bution to the ground state. At finite temperatures we are 
only interested in the second term, which in the scaling 
limit can be written as 



F l 



( sm7 r-^)ln(l + rtf(0) 



2lT T^i J-oo cosh (C - In f ) + costt-^2 



(3.33) 



The evaluation of F l involves the functions rf^ (see Fig. 
|]). The different scenarios possible are very similar to 
those of the multichannel Kondo model. As long as 
fi < /, the impurity remains overscreened and the tem- 
perature exponents are the same as for the \i = 1 case. 
In this case, the residual entropy is 



in- 



nt 



sin 



f+N 



f+N 



nti 1 «n *' 



(3.34) 



(Notice that when N = 2 this reduces to the multichannel 

s j n tt(2S + 1) 

result In — . f t 2 ) . Furthermore, it can be easily shown 

sm 7+2 

that SI =5}_ M . 

We have plotted S l in Fig. |, for fixed N + f = 22, 
and several values of 7 = f /N and [i. Only the region 
/i < / is physical in the figure, since these are results for 
the overscreened case. We can see that for fixed (even) /, 
the largest value of the entropy corresponds to fi = f /2. 
Also, for fixed n, the entropy is the largest around 7 ~ 1, 
and decreases as 7 moves away from 1. 




FIG. 5. Overscreening residual entropy, Sl < t, for an im- 
purity in a totally symmetric representation, \i, of SU(N), 
for N + f = 22. The lower triangular region (in black) is 
unphysical 

If \x = / the impurity becomes completely screened: 
X* and G % /T become constant, and there is no residual 
entropy. Finally, if fi > /, the impurity is underscreened. 
The dominant contribution to the free energy from the 
spin sector is of the form 



F' 



-Tin- 



n^/^sinhrf 



Ilr=i sinnr T Ilr=i sinhr|r 
When H = 0, the residual entropy is 



S l = ln 



(»-f + N-l)\ 
( M -/)!(/V-l)! 



In 



H-f+N-1 

N - 1 



(For N = 2, the residual entropy is In ( ^ ^ ^ 
ln(/z -/ + !)). 



IV. NUMERICAL ANALYSIS 



7+lv 



1G 



A. Procedure 



1. Entropy 



We have solved the TBA equations by iteration, using 
a procedure inspired by the work of Rajancd. For the 
levels /, which do not have a driving term (see Fig. 
H), we have used (^2,^3), as the starting point, since 
it is more convenient to use the kernel G(£). We have 
dealt with the equations that have a driving term by 
introducing two sets of auxiliary functionsa, 

h r (0 = Q(ln(l + tf+i) + + »#_!)), 



so that 



In 77} = 



sin 



(f)+M- + ^) + Q r - 



We have introduced a cutoff A in the integrals involved, 
we have taken ln(l + r^) to be constant for |£| > A, and 
evaluated the integrals in those intervals analytically. For 
|£| < A we have replaced the integral with a sum using a 
Gaussian quadrature ruleo. 

The results that we present in this work correspond to 
zero magnetic field, which means that xq = H/T = 0, 
and the functions rf n depend on £ only. Thus, the task of 
obtaining thermodynamic properties is greatly simplified. 
First of all, the impurity contribution is given by (3.6) 



— Y 

9=1 



(sin^)ln(l + r#(0) 



cosh (£ 



+ COS 7T 



N-q ' 

N 



with rf n independent of T . The entropy and specific heat 
are obtained by taking derivatives of F l with respect to 
the temperature, which can be done analytically when 
xq = 0, and then performing the integration numerically. 

In order to calculate the susceptibility at zero magnetic 
field, x% we derived a second set of TBA equations for 
the functions 



K(0 



dx 2 Q 



2:0=0 



following DegrangesEH. This system is solved as the pre- 
vious one, and the magnetic susceptibility is given by 



X |xo=0 



—Y 

q=l 



d 2 F l 
dxl 



cosh (£_ - In 2a) + cos n^p- ' 



B. Results 



We start by discussing the impurity contribution to the 
entropy, S l . In Fig. ^ we have plotted S l as a function of 
T, for different values of N, /, and impurity spin fi. The 
horizontal axis is on a logarithmic scale. The vertical axes 
have different scales for the different N. The first thing 
to notice is the crossover around T ~ To. For T ^> Tq, 
S % is that of a free spin characterized by (x and N . When 
H = 



S l = In 
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FIG. 6. Impurity contribution to the Entropy as a function 
of T, for different values of N, f, and impurity spin, [i. Notice 
that different scales are used for the different N. 

Below the crossover region one can see the quenching 
of the degrees of freedom due to the interaction in the 
decrease of the value of the entropy. The qualitative be- 
havior of S l in the region T -C Tq depends only on the 
relation between fi and /. When / = ji, there is com- 
plete screening, and S l = 0, as can be seen in the curves 
fi = f = 1, fi = f = 4. For fi > /, the impurity is 
not completely screened, and there is effectively a left- 
over free spin fi — f, as can be seen in Fig. 6 for / = 1. 
There, the T <C To entropy for fi < f corresponds to 
the T > To entropy for \i — 1. Finally, when fi < /, 
overscreening takes place: even though there are enough 
electrons to form a singlet with the impurity, the low 
temperature behavior is characterized by an object with 
complex internal structure, and an anomalous S l . Such 
behavior can be seen in the curves fi — - 2, fi — 3 for / = 4, 
and in all the curves for / = 8. Notice that for fi = 1 
and fi = / — 1, the curves converge to the same value, 
as we had already seen in the asymptotic analysis. Fur- 
thermore, the overscreened fixed point has an anomalous 
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residual entropy irrespective of the value of N, indicating 
its Non-Fermi liquid nature. 

A more detailed picture of the behavior of S 1 for N = 4 
is displayed in Fig. [?]. It is worth noticing that in the 
underscreened cases, the effective spin is fJ, — /. Also, 
there might be situations where the residual entropy of 
the overscreened case is larger than that of the under- 
screened case. Such is the case for / = 5, fx = 3,4, as 
compared to fi = 6. 
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Next, we compute the contribution to the specific heat. 
Results for different values of the parameters are shown 
in Fig. |[ The largest maximum of C v corresponds to 
fj, = f. Also, the size of the curve grows with N. 

We have also evaluated the subleading contribution to 
the linear coefficient of the specific heat, j l = C v /T for 
f = 2fj, = 2, N > f, and we have plotted it in Fig. 
|9|. The points fit power-law curves with exponents (N — 
f)/(N + f), derived previously (see Eqs. ( |3~29| ), ( pp| ), 
and (3.32)). . This is another clear indication that for 



N > f, the overscreened cases are not Fermi liquid fixed 
points. 



f=2 




FIG. 7. vs. T for JV = 4, / = 1, 3, 4, 5, and fi = 1, . 



2. Specific Heat 
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FIG. 8. Impurity contribution to the specific heat, C v , for 
different values of N, f, and /j, 



FIG. 9. Subleading contribution to j(T) = C v /T, for 
/ = 2, fi = 1, and N > f vs. T, in a Log-Log graph. The 
symbols correspond to the numerical calculation. The lines 
correspond to power- law fits with exponents (N — 2)/ (N + 2) . 



5. Magnetic Susceptibility 

Next, we have studied T^ 1 for different values of the 
parameters, and plotted the results on Fig. As with 
the entropy, the qualitative behavior depends only on the 
values of /i and /. The difference in behavior between un- 
derscreened and overscreened cases becomes more clear 
here: whereas the magnetic moment is partially quenched 
in the former case, the overscreened case is characterized 
by a totally quenched moment, even though there is a 
non-zero residual entropy. This can be seen in Fig. ^ 
for the curves with / = 1 and / = 4. 

The magnetic susceptibility is plotted in Fig. |ll[ The 
curves with / = fj, have a constant % l at low-T and for 
N > 2 they have a maxima near T ~ To- This is a 
special feature of the completely screened case. We see 
that in the overscreened case with N > /, the suscepti- 
bility tends to a finite value as T — > 0, while it diverges 
when N < f. When N = f, the divergence is loga- 
rithmic, whereas it is power-law for N < f, with an 
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exponent /3 > — 1. The largest divergence corresponds 
to the underscreened case with 1/T behavior. All these 
result s coin ci de w ith tho se of the previous analytic study, 
Eqs. (HH) , jOcI ), and (§J|). 









1=1 


f=2 






^' 


_ — n=l 
-- fl=2 
H=3 







- m /: ;....., 




/ ; 












y*' 






~~ .X 




/'... 





/ 

J/ 1 


/' 

/ ,. 






// 




— / 













\ V 




H=f=l 




I I iT|*T-Tl 




\ i— ► 


N=2 


H=f=4 


\\\\ 


N=4 

N=6 

N=8 



10 
T/T„ 



FIG. 12. x i (r)/x*(0) vs. T for different values of N, in 
two completely screened cases, fj, = f = 1, and /x = / = 4. 
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FIG. 10. vs. T, for different values of /, N, and fi 
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FIG. 11. Impurity contribution to the magnetic suscepti- 
bility, x 1 vs. T, for different values of /, N, and fi. 

In Fig. |l2] we show \ % f° r different values of N for the 
cases /i = / = 1 and // = / = 4. We have rescaled the 
curves dividing by x*(0)- It is quite apparent that the 
behavior is the same in both cases. 




FIG. 13. x i vs. T for AT = 4, fj, = 1, and / = 1,2, 4, 6, 8. 

In Fig. |l4] we consider N = 8, f = 7 and p = 1, 7. 
Even though x*(0) is finite, it is clear that the behavior 
of the overscreened case is quite different from that of the 
completely screened case, and that the subleading terms 
have an important contribution below the crossover tem- 
perature. 

We can also see the power-law behavior of the sub- 
leading term of x 1 m Fig. 15. The values for the ex- 
ponents agree with the values obtained analytically, i.e. 
(N-f)/(N + f). 
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N=8, f=7 




FIG. 14. x 1 vs. T for N = 8, / = 7, and ^ = 1, ...7. 
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FIG. 15. Subleading T dependence of x* for N = 8, 
f = 2, 7, /i = 1. The lines correspond to power law fits. 



4- Wilson Ratio 

We have calculated the Wilson ratio, defined as 

TX 1 



R 



- 1)m 2 b Cj 



(4.1) 



The quantity R has a well-defined meaning only for T = 
0. However, in Fig. [if] we have plotted the quantity 
R{T), to show the difference between the N < f and 
the N > f sectors. In the former case, (N = 2, / > 
1), the value for the overscreened case is much larger 
than the value of the completely screened case (notice 



the difference in vertical scales) , whereas in the latter case 
(N = 8), the curves converge to the completely screened 
value. Notice that for N = 4 there is a change in behavior 
as we go from / < N to / > N. 
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FIG. 16. R(T) as a function of temperature for N = 2,4, 8, 
M= 1, -,/ an d / = 1,3,5. 

Next we plot the values of R for completely screened 
cases, [i = f (Fig. [I?]). We see that the values obtained 
fit the function 



R = 



N(N + f) 
N 2 -l ' 



(4.2) 
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FIG. 17. Wilson Ratio R for /i = /. The points have been 
obtained from the numerical solution. The lines correspond 
to fits with the function N(N + f)/{N 2 - 1). 

We have also obtained values of R for the overscreened 



20 



case (/i < /), in Fig. ffq. There are clear differences 
between the / < N and the f > N cases, as we have 
already pointed out. For / < N, the valu e of R coincides 
with the value for fi — /, and agree with (4.2). The dom- 
inant contribution to R comes from the constant terms 
in x l ancl Cy/T. For / > N, R contains mainly the 
coefficients of the divergent parts, and have a different 
functional behavior. 



(when fi < 6), characterized by a SU(6) flavor symmetry. 
The pairs of curves (fi, 6 — fi) merge. Notice also that 
fi = 6 is completely screened. There is a second crossover 
around T ~ T4 to a region characterized by SU (4) flavor 
symmetry for fi < 4 and by SU(2) flavor symmetry for 
4 < fi < 6. Only the curves for fi = 1 and fi = 4 — 1 = 3 
coincide now. Also, the value of the residual entropy for 
fi = 1, is that of an effective fi = 1 in a SU(3) x SU(2) 
model. The entropy for fi = 4 goes to with T since the 
system becomes screened for T <C T4. 
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FIG. 18. Wilson Ratio 7? for fi < f. See the caption in the 
previous figure. 

5. Channel Anisotropy 

We end the discussion of the numerical results by show- 
ing an example of channel anisotropy. We have taken the 
case N = 3, with original flavor symmetry SU(6), broken 
down to ££7(4) x SU{2) (Fig. |9|) Two scales appear in 
this problem: Tq, and T4. Accordingly, the TBA equa- 
tions have driving terms at n = 4, 6. 




FIG. 19. Diagrammatic representation of the TBA equa- 
tions for a system with N = 3, and flavor symmetry SU(6) 
broken down to SU(4) x 517(2). Two driving terms appear 
at n = 4, 6 

We have chosen very small anisotropy, T^/Tq = 10~ 6 . 
In Fig. |o| we have plotted the entropy for different values 
of fi. There are three different regions: When T » 
Tq, the impurity behaves like a free moment. Around 
T ~ Tq there is a crossover to an overscreened region 
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FIG. 20. Entropy vs. T, for fi = 1, 6 
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FIG. 21. Cy/T vs. T, for fi = 1,4,5,6. 

Finally, we have plotted Cy/T for several fi in Fig. 
When T 4 « T << Tq, the behavior of the system 
is characterized by / = 6: increasing value of Cy/T for 
fi < 6, constant behavior for the screened case, fi = 6. 
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Near T4 the curves for \i < 6 have a similar behavior as 
those for free moments (however, such behavior is not 
found in the curves for the magnetic susceptibility). We 
can see three different behaviors in the region T C T4. 
For [i — 1 the curve diverges with a power law, as one 
would expect for f e f j = 4 > N = 3. The curve for fi = 4 
is flat, since the impurity is screened. Finally, for ji = 5, 
the curve increases slowly, converging to a constant value, 
as one would expect from an effective f e f f = 2 < N = 3. 



V. PHYSICAL REALIZATIONS OF THE 

SU(N) x SU(F) MODEL 

In the present section we shall concentrate on the 
possible physical realizations and applications of the 
SU(N) x SU(f) Coqblin-Schrieffer model. 



In Subsections. V A, V_B, and V C we analyze the so 



called N-level system (NLS)-. model, a generalization of 
the two-level system modelJ3 by means of a systematic 
1/f expansion. In this context the flavor degeneracy is 
associated with the physical spin of the electron. First 
we establish a mapping of the NLS model to the multi- 
channel Coqblin-Schrieffer model (MCCS model) by ana- 
lyzing the low-energy fixed point of its scaling equations. 
While our procedure gives a systematic expansion only 
for the case / > N we shall argue that the same mapping 
should apply for the cases f < N. 

In the limit / > N we are able to determine the full 
operator content of the fixed point. This enables us to 
calculate the scaling of the differe nt ph ysical quantities 
at low temperatures in Subsection V C . As we shall see 



there are some subtle differences between the two models 
and while most of the physical quantities show the same 
dependence, the scaling of the specific heat may be dif- 
ferent. The origin of these differences will be discussed 
in detail. 



Finally, based on the res ults of Subsections V A - VC , 
we discuss in Subsection VD some physical systems 
providing possible candidates for the realization of the 
MCCS-model. 



A. The N-level system model and its low-energy 
fixed point 

The NLS model has been constructed as a generaliza- 
tion of the two- level system modelBo to describe the tun- 
neling of a heavy particle among TV not necessarily equiv- 
alent positions labeled by a — {1,..,N}, and strongly 
coupled to the conduction electrons. At low tempera- 
tures the motion of the heavy particle can be described 
by the effective Hamiltonian 



N 

0,6=1 



(5.1) 



where xt creates a pseudofermionEa corresponding to 
the heavy particle site a and A afc is the tunneling am- 
plitude between positions a and b. If no external stress is 
present then the diagonal part of A ab vanishes: A aa = 
(a = 1, ..,N), when the N positions are equivalent due 
to the symmetry of the NLS. The electronic part of the 
Hamiltonian and the coupling of the heavy particle to 
the conduction electrons take the general form: 



e c 



- 



enm 

H e — hp = ^ c enmXa^nn'Xb c e'n'm ; 
a,b 7 n,n 



(5.2) 



where the operators cf nm create conduction electrons 
with energy e, orbital quantum numbertl (~ angular mo- 
mentum) n = 1,2,.., 00, and spin m. For the sake of sim- 
plicity the electronic density of states g(e) is assumed to 
be constant go, between the high- and low-energy cutoffs, 
D and —D, independently of the flavor and orbital quan- 
tum numbers. While in the physical case only m = ± is 
possible corresponding to the two different spin direc- 
tions, for technical reasons in the following we assume 
that the electron spin m can take / different values: 
m= 1,..,/. 

This model has a structure similar to Eq. (2.1) but 
there are some important differenc es. The 'spin index' 
a' = 1, ..,N of the impurity in Eq. (2.1) is now replaced 
by the 'site index' a = 1,..,N of the heavy particle. 
Moreover, in the NLS case the orbital index n (replacing 
the 'spin index' a = l,..,N of the conduction electrons 
in the SU (JV) x SU(f) model) now ranges from one to 
infinity since the conduction electrons may have any or- 
bital momentum. Furthermore, the couplings are highly 
anisotropic in orbital indices and no SU(N) symmetry is 
present at this level. Finally, in the NLS model the scat- 
tering is diagonal in the real spin inde x m which plays 
now the same role as the flavor in Eq. (2.1 ). 

The diagonal couplings, V^, describe simple poten- 
tial scattering of the conduction electrons by the heavy 
particle sitting in position a. On the other hand, the off- 
diagonal matrix elements, V^, with a ^ b correspond 
to the so called "assisted tunneling" processes. Here the 
heavy particle is tunneling from one site to another while 
a conduction electron is scattered by it. The combina- 
tion of these two processes leads ultimately to the gener- 
ation of an orbital Kondo effect and a strongly correlated 
ground state.ErB 

In the following we shall carry out a large / analysis to 
determine the low-energy fixed point of the NLS model. 
While our procedure is strictly valid only in the case / > 
N, in the end of the subsection we shall argue that our 
results are very general and they should apply even for 
the cases / < N. 

To carry out a 1// analysis of the NLS model, as a next 
step, we construct the next to leading logarithmic scaling 
equations using a generalized multiplicative renormaliza- 
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tion group technique.^ As discussed in Refs. |l|, |l7| , |l8| , p6| , 
the leading logarithmic equations give the leading term in 
a systematic 1/ / expansion and become exact in the / — + 
oo limit. In the multiplicative renormalization group 
method one exploits the existence of a non-trivial trans- 
formation in the space of the Hamiltonians, D — > D/, 

Kn> -» y'nn', and Aab A ' ab th at leaves the pseud- 
ofermion Green's function Q ab and the pseudofermion- 
conduction electron vertex function TtL,(w) invariant: 



G{uj, V, A', D') = A Q{uj, V, A, D) A+ , 
IV, V, A', D') = [A + ]- 1 T(u;, V, A, D) A' 1 



(5.3) 



In these equations A denotes an N x N matrix inde- 
pendent of the energy variables, u> and T, and acting in 
the site indices: A = A a HV',A',D'/D). By means of 
the transformation Eq. ( [5.3| ) one can generate effective 
Hamiltonians that describe the system's behavior below 
the energy scale D'. The generated effective Hamiltoni- 
ans usually (in a renormalizable theory) turn out to be 
much simpler than the original one. 





a. 




b. 



FIG. 22. The leading logarithmic vertex and self-energy 
diagrams generating the next to leading logarithmic scaling 
equations. Continuous and dashed lines represent the con- 
duction electron and pseudofermion Green's functions. 



To make use of the invariance property Eq. ( p.3[ ) one 
first has to construct the lowest jorder vertex and pseud- 
ofermion self-energy correctionsa arising from the dia- 
grams in Fig. 22. a and b, 



-/In — (S ab u tr{« c V c } ~ tr{i; ac A c V fc }) 
D 



QoL ab = H ah - In — ([v ac , v cb ] - f tr{v ac v db }v cd ^j , (5.4) 

where go is the density of states at the Fermi level, and 
a matrix notation has been introduced for the dimen- 
sionless couplings go V^ b , — > v ab . The symbol [ , 
stands for the commutator, the trace operator tr{...} is 
acting in the electronic indices, and a summation must 
be c arrie d out ove r re peated indices. Then, substituting 
Eq. ( 5A ) into Eq. ( |5.3|) and reducing the bandwidth D by 



an infinitesimal amount one can deduce the infinitesimal 
renormalization group transformations for the couplings: 



dA a ' 



= - l -f 
dx 2 J 



tr{v ac v cd }A db + A ac tr{v cd v db } 



-2tr{w ac A c V } 



dv ab 
dx 



k ac ,v cb ] + ^f(2tr{v ac v db }v 



cd 



j_„ r ac crfi db aci_ r cd dbi \ 

tr{v v \y_ — v tr{v_ v }J, 



cd„.dbi 



(5.5) 



(5.6) 



where the dimensionless scaling variable x = ln(Do/D) 
has been introduced, Do being the initial (real) band- 
width cutoff of the model. 

These scaling equations have to be solved with the 
boundary condition that the couplings are equal to their 
bare values at x = 0, and they loose their validity if the 
reduced bandwidth D becomes smaller than any small- 
energy scale present, T, u>, A. Note that, up to the next 
to leading log arith mic order, the splittings A ab do not 
occur in Eq. (5.6) explicitly, and they provide only a 
low-energy cutoff for the scaling. To be explicit, there 
is an energy scale, D* = T* that we call the freezing 
temperature, where the renormalized splitting becomes of 
the same order of magnitude as the reduced bandwidth: 
A ab (D*) ~ D*. Below this energy scale the orbital mo- 
tion of the NLS is usually frozen out (see the discussion 
in the end of this Section), and the couplings may be 
replaced by their values at T* . 

For the moment let us forget about Eq. (|5. 5|) a nd con- 
centrate on the scaling of the v ab, s, Eq. ( |5.6| ). This 
equation cannot be solved generally, but one can con- 
vince oneself very easily that if the 'assisted tunneling' 
matrix elements v ab (a ^ b) do not vanish then the 
electron-NLS ffmiplings start to increase and lead to a 
Kondo effectfl'ta The scaling of the norm of the cou- 



plings, V , 



is shown in Fig. E23| for a symmetrical 



six-state system, where the coupling constants have been 
estimated using similar methods as in Ref. ||. As one 
can see, a Kondo effect occurs around the Kondo scale 
T K ~ Do e~ x - - WK, where x c = In (D /T K ) denotes 
the value of the scaling parameter at which the crossover 
from weak to strong coupling occurs. Our numerical in- 
vestigations for various model parameters and different 
values of N show that the structure of the stable low- 
temperature fixed point the couplings scale to is inde- 
pendent of the initial couplings and only depends on the 
value of N as long as no special symmetry is assumed for 
the v ab, s. 

In what follows we shall show that this stable low- 



energy fixed point of Eq. (5.6) has the structure of the 
defining representation of the SU(N) Lie algebra. To 
be precise we first observe that the operators O a ~ 
5 ab ^2 c v cc are invariant under scaling. Therefore the 
v ab, s can be divided into two parts, v ab and M ab where 
So = ano - M. ab is built up from the previously men- 
tioned constants of motions, O a . Then as we shall see, 
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at the stable fixed points of Eq. ( |5.6| ) the v ab, s can be 
written as 



(£ a6 )fp = 7 ( "r 



1 /-L ah 



(5.7) 



/ V o o, 

where the L ab, s satisfy the standard SU(N) Lie algebra, 
[L ab , L cd ] = 5 ad L cb - 5 cb L ad . (5.8) 



and are unitary equivalent to the defining representation: 
L nn> ~ S n' S n -jS ab &nn'- This statement is also demon- 
strated in Fig. E3] where the scaling of the 'algebra coef- 



£cb~ad\ 



IS 



ficient' a = Z a ,b,c,d 11/^,^1 - fS ad v c " + f5 cb V 
shown, measuring how well the fixed point algebra ( |5.8| ) 
is satisfied. As one can see in Fig. ^3] for D <C Tk the 
algebra coefficient a vanishes and therefore, in an ap- 
propriate basis, the v ab, s really simplify to the form in 



Eq. (5.7) 



3 

B 




5.00 10.00 15.00 

x = In ( D /D ) 



FIG. 23. Scaling of the norm of the dimensionless cou- 
plings, u — ^2 \ \v ab \\ (dashed line), and of the algebra coeffi- 
cient a (continuous line) for a 6-state system with / = 2. 

E q. (|5.7|) means that at the fixed point v ab is given by 
Eq. ( |5.7| ) and apart from some potential scattering term 
the fixed point effective interaction can be written as 



H c s — Vo 



E 

a . 6 



Xt c e bn 



iXb 



(5.9) 



which is the same as the interaction term in Eq. (2.1). 
Note that while the initial model was very asymmetrical 
in the orbital space, at the fixed point only N conduc- 
tion electron angular momentum channels are coupled to 
the NLS and the fixed point effective Hamiltonian shows 
already an additional SU(N) symmetry in the orbital 
sector as well. These statements are not true away from 
the fixed point where various kinds of irrelevant oper- 
ators couple the NLS to the electrons and coupling to 
the other orbital channels is also relevant. The effective 
Hamiltonian is completely symmetrical in the NLS site 



index, which means that, e.g., the amplitude of assisted 
tunneling from site 1 to 6 in Fig. ^ is the same as the 
nearest neighbor assisted tunneling amplitude from site 
1 to 2, despite of their different geometrical position. 

As it will become obvious in the next subsection, the 
analysis above is based on the possibility of a system- 
atic N I f expansion. Therefore, it is strictly valid in the 
/ > N case. However, one has several arguments that 
the effective Hamiltonian Eq. (5.9) is also adequate for 
the N < f cases. First of all, in the case N = 2 cor- 
responding to the simpler case of multichannel Kondo 
model, it is well-known that for / = 1 and / = 2 the 
spin anisotra»¥ of the couplings is irrelevant around the 
fixed point pc3 which has the same SU(2) structure as 
Eq. (gj). Furthermore, for / = I but arbitrary N one 
can easily prove following similar lines as Nozie res and 
Blandiro that the isotropic fixed point, Eq. (5.S), is sta- 
ble against spin (orbital) anisotropy. These observations 
together with our results for the / > N case make it 
highly improbable that for 2 < / < N the NLS model 
would have a stable fixed point different from the one 
discussed above. 



B. Stability analysis of the SU(N) x SU(f) fixed point 
of the NLS model in the large / limit 



The statement that v ab is a fixed point of Eq. (5.6) 
is trivial. However, we also want to prove the stability 
of this fixed point analytically and find the irrelevant 
operators (determining the low-energy behavior of the 
model) around it. To this end we write the deviations 
from the fixed point in the form: 



Sv ab = 



(t 



n ab fab 
ba\-\- ^ab 



(5.10) 



where the couplings g ab , t ab , and [i ab are N x N, N x oo, 
and oo x oo matric es, respectively. Substituting this 
expression into Eq. (|5.6| ) one obtains the following lin- 
earized decoupled scaling equations: 



dft 



ab 



d.r 
dg ab 



f 



( S ab^dd _ N ab} 



(5.11) 



1 



, =~-A\L ad ,Q db } + [Q ad ,L db ] 
dx J 



_Lf 2 5 ab Q dd + 2L cd tr{g ac L db + L ac g db } 

Z / v 



2/ 

2N g ab - L ac tr{g cd L db + L cd g db } 
tr {g ac L cd + L ac g cd }L db } 



dt ab 
dx 



f 



jjidfdb 



f 



(5.12) 
(S ab t dd - Nt ab ) . 



(5.13) 

The solution of Eq. ( |5.11| ) is trivial, since the operator /i ab 
can be decomposed as [i ab = (fi ab -5 ab ±n cc ) + (6 ab j r fi~ c ), 
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where the first operator scales like 



- x A 



with a di- 



mension A = A s i = N/ f , while the second is marginal 
with A = 0. The detailed analysis of the other two equa- 
tions is much more complicated, but one can still find 
their exact solutions due to the simple structure of the 



L s Here we only briefly discuss the results of this anal- 
ysis. 



It turns out that Eqs. ( 5.12 ) and ( 5.15 ) have an infi- 
nite number of zero modes that can be divided into two 
distinct classes. The first type corresponds to potential 
scattering off the NLS and can be written as 



S ab Sv 



(5.14) 



where Sv_ denotes an arbitrary oo x oo Hermitian matrix. 
The rest of the zero-modes can be identified with the 



generators Sv® cn of the unit ary transformations of the 
SU(N) Lie algebra, Eq. (§H 



^gon 



(t ba ) 



gen 




(5.15) 



More precisely, the generators g^ n and t^ n can be shown 
to satisfy in first order the equations: 



r abj-cd j 
±L igen ~~ Li 



cd±ab 



cadj-cb 
Cp-en 



ccb^ad 
igcn i 



\L ab , g cd 1 + \L cd , g ab 1 = S ad p cb 

— — gcn J L — — gcn J —gen 



5 cb t 



(5.16) 



; ab 



from which it follows that the operators L 
Sv °g n satisfy the same Lie algebra as the original L 



jab 

ab }„ 



[L ab ,L cd ] 



S ad L cb 



5 cb L aa 



(5.17) 



The o 's turn out to be the generators of the uni- 

— gen ° 

tary transformations in the iV-dimensional electronic 
subspace where the Lie-algebra Eq. (5.8) is realized, 
while the t ab, s correspond to the rotations of this N- 
dimensional subspace. 

All the other eigenoperators around the fixed point 
can be shown to be irrelevant. Very surprisingly, at least 
in the large / limit, the leading irrelevant operators are 
quite different from the leading irrelevant operator of the 
SU(N) x SU(f) model d5~l), both in their structure and 



in their scaling dimension. They are living in the sector 
t ab and they can be written as 



c ab 
OV i 





(C ba )+ 



c ab 





(5.18) 



where the C^'s satisfy J2 a C Z = and J2b( C 2™ ~ 
Cf m ) = with a, 6, c = 1, ..,7V and n = 2V+l,iV+2, ..,oo. 
These operators have a dimension 



A, = 



N -1 



f \f 2 



(5.19) 



We remark at this point that the operators (5.18) do 
not exist in the two-level-system model, which is there- 
fore completely equivalent to the corresponding SU (2) x 



SU(f) modelB As we shall see, these operators do not 
give a contribution to physical quantities like the resis- 
tivity or the impurity susceptibility but they influence 
the thermodynamic behavior of the model. We stress at 
this point, that their existence is strictly proven in the 
/ — > oo limit. They are very probably present even in 
the N < f case but it is an open question if they survive 
in the N < f limit. 

The impurity resistivity will be shown to be dominated 
by the subleading operators 



c ab 



Q 



ab 



o 

o s ab 



(5.20) 



where the matrices Q ab and S ab satisfy Q aa = S aa = 
and Qfc = Q^ b - These operators have a dimension 



N 



A s , = — + ^ — 



N 2 



P 



(5.21) 



and the operator ( |5.9| ) considered in the SU{N) x SU(f) 
model is also one of them. Furthermore, one has other 
even more irrelevant operators in the t sector of Sv_ with 
a dimension A ss i = (N +1)// which give a subleading 
contribution to the physical quantities calculated. 

In the previous considerations we did not take into ac- 
count the presence of the splitting A ab of the NLS. As we 
discussed already, this splitting results in the appearance 
of another low-energy scale, T*. Below this the NLS can- 
not jump freely between its N different positions. Since 
usually the ground state of the NLS is non-degenerate in 
most cases a Fermi liquid state develops. In other words, 
the non-Fermi liquid SU(N) x SU(f) fixed point is un- 
stable with respect to the splitting that usually drives 
the system towards a Fermi liquid-ffround state. 

It has been argued very recentlycSl that in special cases, 
due to some dynamical Jahn- Teller effect, e.g., the hop- 
ping amplitude A ab might pick up and additional Berry 
phase, which could then result in a degenerate ground 
state with degeneracy N' . Then the effective Hamil- 
tonian at very low temperatures would be, of course, 
an SU(N') x SU(f) exchange model, and in the region 
T <^T* all our previous considerations hold with the re- 
placement of N by N'. Unfortunately, this Berry phase 
scenario is not very probable and therefore we expect the 
non- Fermi liquid behavior can be only observed in the re- 
stricted temperature (energy) range T* < max{T, w} < 
Tk, i.e., when the freezing temperature is small enough. 

Therefore it is very important to determine the real- 
istic values of the freezing temperature. We estimated 
the fr eezin g te mper ature by solving the scaling equa- 
tions (5.5) and (5J3) numerically for the same symmet- 
rical 6-level system as in Fig. In this case the di- 
agonal matrix elements A aa vanish by symmetry. As 
one can see from Fig. |24|, for a realistic NLS the renor- 
malization of the hoppings A ab is huge, and the situa- 
tion T* < max{T, lo} < Tk can be reached quite eas- 
ily. We note at this point that in our Hamiltonian we 
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also neglected the contribution of two-electron scattering 
around the fixed point, which might be also reLevant in 
the immediate neighborhood of the fixed point £3 How- 
ever, these have a very small amplitude and they are 
scaled downwards in the first part of the scaling, D > Tk- 
Therefore, most likely their effect can be neglected com- 
pared to that of the splitting A ab , which provides the 
dominant—mechanism driving the system to a Fermi liq- 
uid state.LJ 




0.00 5.00 10.00 15.00 20.00 

x = In ( D /D ) 

FIG. 24. Scaling of the dimensionless hopping amplitude, 
A 12 /Dq for the same 6-state system as in Fig. Inset: 
Numbering of the sites of the 6-state system. 



C. Scaling of the physical quantities of the NLS 
model in the large / limit 

Now we turn to the calculation of the physical quan- 
tities. In this Subsection we shall determine different 
thermodynamic quantities and the conduction electrons' 
scattering rate, 1/r, which is directly proportional to 
the impurity contribution to the electrical resistivity, 

Rimp(T). 









d. 



e. 



FIG. 25. Diagrams generating the l// 2 corrections to the 
free energy. The crosses denote the counterterm. 

To calculate a general physical quantity one should also 
calculate the renormalization coefficient A in Eq. (5.3), a 
non-trivial task away from the fixed point. However, one 



can easily convince oneself that the factors A and A^ 1 in 
the free energy corrections in Fig. and the electronic 
self-energy corrections in Fig. g2|c cancel exactly, and 
therefore these are scale invariant, and ca n be calcu late d 
by considering only the scaling equations ( |5.5| ) and (|5.6|). 

For the sake of simplicity let us assume first that the 
highest low-energy scale is given by the temperature. To 
calculate a physical quantity at a temperature T we ap- 
ply a renormalization group transformation (5.3) with 
D = Do and D' = T . Then in the new Hamiltonian all 
the logarithmic terms vanish since \n{D' /T) = 0, and the 
different physical quantities are exclusively given by the 
non-logarithmic contributions of the corresponding dia- 
grams. For a scale invariant quantity like the free energy, 
e.g., this implies that 



F imp (D Q ,T,v% b 1 A a a b ) = 
F imp (T,T,v ab (ln^),A ab (\n^)) 



(5.22) 



where on the right-hand-side no logarithmic corrections 
appear, but the renormalized couplings have to be used. 

Therefore, in order to calculate the scaling behavior 
of the thermodynamic quantities our task is to deter- 
mine the non-logarithmic parts of the different free en- 
ergy diagrams. Since the fixed point couplings v® b are 
proportional to 1// up to 1/ f 2 order only the diagrams in 
Figs. ^5].a-d. contribute. However, these diagrams con- 
tain divergent contributions originating from the finite 
part of the self-energy diagram in Fig. p2lb. These spuri- 
ous divergences can, be handled by a standard renormal- 
ization procedure,E3 by adding the following counterterm 
to the Hamiltonian: 

Hcount = / 2D In 2 xtXc tr{v a b v bc ] . (5.23) 

This counterterm can also be interpreted as a renormal- 
ization of the ba re p arameters of the model, which should 
be used in Eq. (5.3) as the initial conditions. Then the 
counterterm contributions in Fig. |2^.e cancel all the spu- 
rious divergences, and after a tedious calculation one ob- 
tains for the non-logarithmic part of the free energy: 



imp 
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In AT- 



37V 



(tr{v ab v bc v ca - v ca v bc v ab }) 



2^2 



(tr-jV 



ab n cd 



bc„da 



<— } l r{r ' i 



2N 

tr{v ab V bc }tr{v cd V da }) 



} 



(5.24) 



Note that diagram (a) in Fig. ^5] is proportional to T 2 /D, 
and it does not give a contribution in the scaling limit. 

Su bstit uting the fixed point couplings Eq. (5.7) into 
Eq. ( 5.24 ) the fixed point entropy can be calculated as 



Sin 



dF h 



dT 



In AT 



N 2 



1 7T 



P 



G 



(5.25) 



which is just the expanded version of Eq. (3.18). Note 
that Eq. (5.25) gives the NLS contribution to the entropy 



2G 



only in the region T* < T < T K . Below T* the motion 
of the NLS is usually frozen out and the impurity entropy 
tends to zero, as expected for a Fermi liquid state. 

The scaling of the free energy in the region T' < T « 
Tk- can be determined by expa nding the v ab, s around 
their fixed point values like Eq. ( 5.1 C| ) and substituting 
them into Eq. (5.24). It turns out that similarly to the 
multichannel Kondo and the two-level system caseE3'E£l 
only the second order terms in Sv ab contribute, and there- 
fore in the temperature range T* « T < Tjc in leading 
order the free energy and the specific heat scale as 



-Pimp ~ T | ^ 



T 

K 



2A, 



2(jV — 1) 

Tk~ 



(5.26) 
(5.27) 



Below T* the Free energy generally shows a Fermi liq- 
uid behavior. This scaling behavior does not agree with 
the one obtained in the BA solution of the exchange 
model. However, we have to remark at this point, that 
according to our estimations the amplitude of the sub- 
leading operators in Sv ab is larger than that of the lead- 
ing irrelevant operators. Therefore, one expects that 
there is a substantial energy region where the sublead- 
ing operators dominate, and eventually it is also possible 
that they dominate the scaling of the free energy in the 
whol e reg ion T* <C T <C Tk- Then the exponent Ai in 
Eq. ( |5.26 ) should be replaced by A s i and one obtains a 
scaling Ci mp ~ T 2N ^ which is in 1/f order completely 
identical with the Bethe ansatz and conformal field the- 
ory results for the SU(N) x SU(f) model. 

One can also easily determine the scaling of the split- 
ting susceptibility xa = d 2 Fi mp / dA 2 at T = for small 
A's, where now A denotes the characteristic value of 
the splittings A ab . Investigation of the free energy di- 
agrams Fig. p5| shows that the 'splitting magnetization' 
Ma = dF im p/dA should be of the form: 



values. Since m(l,v^ b ) is just a constant, the scaling of 

Ma is the same as that of the factor Za (§r,v ab ). For 
very small A's the scaling of can be easily deter min ed 
from the fixed point form of the scaling equation fl5.5|) 



dA 



ab 



N 



-A 



ab 



dx f 
and one obtains in leading order in 1/f: 



(5.31) 



(5.32) 



in agreement with ( p.29|) and the conformal field theory 
resultsO In higher order in 1// one also has to ta ke int o 
account the renormalization of the splitting in Eq. (p. 32 ), 
a* _ aV(i-Asi) and one obtains with A=- - ~ 
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JV 

/ 



N+f 



\ - 

P 



M A ~ A~ ~ A N ' S , (5.33) 

which is the exact result. 

Finally, we discuss the scaling of the electronic scatter- 
ing rate, which we determine from the imaginary part of 
the electronic self-energy in Fig. ^.c. By assuming a fi- 
nite impurity concentration rii and averaging over the po- 
sition of the impurities and the orientation of the incom- 
ing electrons we obtain for the average scattering rate: 



(-) = 2nm(2D ) ^ tv{v ab v ba }. 
t N 



(5.34) 



Note, that the factor Do arises from the inverse density 
of states Qq 1 and is invariant under scaling. Substitut- 
ing into this equation v ab — vf b + Sv ab we see immedi- 
ately that the leading irrelevant operators do not give 
a contribution to the electronic scattering rate which is 
dominated by subleading operators and scales like 



M A =m(-,v ab ) 



(5.28) 



The important point is that A is not scale invariant, but 
it rather behaves as 



A' 



A 



(5.29) 



where the factor should be determined by integrat- 
ing Eq. (5.5). As a consequence, Ma is not scale in- 
variant either, and it has to be rescaled under the RG 
transformation by the factor Za- Therefore, ap plying 
the renormalization group transformation to Eq. ( 5.2§| ) 
with D' = A* = T* we obtain: 



xm(l,vf b ) , 



(5.30) 



where we assumed that A* <C Tk and thus the scaled 
couplings v ab (D') can be replaced by their fixed point 



~ T N /f (w = 0) 
A ' (T = 0) . 



UJ 



(5.35) 



In higher orders this result should be replaced by 1/r 

T N/(f+N) and 1/r ^A7(/+A0. 



D. Discussion of the possible physical realizations of 
the NLS model 

The simplest possible realization of the SU (N) x SU(f) 
model is given by substitutional impurities jji-metals. 
These impurities may form tunneling centersEju which 
then interact with the conduction electrons' band. An ex- 
ample of such a system is given by Pbi_ x Ge x Te!El The 
alloy PbTe is a narrow gap semiconductor, but usually 
because of some intrinsic impurities it becomes metallic 
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at low temperatures. Since the Ge 2+ ions are smaller 
than the Pb 2+ ions and they are also attracted by their 
nearest neighbor Te 2 ~ ions, they form 8-state syst ems, 
and according to the our discussions in Section V A they 
would be good candidates for the SU (8) x SU (2) model. 

However, while an unambiguous logarithmic anomaly, 
has been observed in the resistivity of these materials 
no non-Fermi liquid behavior has been detected. There 
may be several reasons for that. According to the results 
of the BA calculations in the case N > f the non- Fermi- 
liquid corrections are subleading, and the physical quan- 
tities have in leading order a Fermi liquid-like behavior. 
The subleading low-temperature behavior of these alloys 
has never been analyzed and the original measurements 
do not seem to be accurate enough to extract such de- 
tailed information. We are also not aware of any mea- 
surements of other physical quantities like the specific 
heat in the interesting concentration domain. Further- 
more, PbTe has very complicated properties: it has a 
soft phonon mode that drives the system trough a fer- 
roelectric phase transition as a function of the Ge con- 
centration, and strong spin-orbit scattering which most 
likely spoils the SU(2) symmetry of the electron spins as 
well. Moreover, the measurements have been carried out 
at relatively large Ge concentrations, where the interac- 
tion of the NLS's can not be neglected any longer. 

It seems that in order to observe the non-Fermi liq- 
uid scaling much more accurate measurements should be 
carried out at even lower temperatures and lower Ge con- 
centrations for several physical quantities and lower Ge 
concentrations. One could also try to find a better can- 
didate. Since in the case of PbTe the formation of the 
NLS's is induced by the ionic attractions, we think that 
experimentalists should search among multicomponent 
metals, where interstitials can occur in the material. 



SU(f) model could be possible by means of nanotech- 
nologies. In Fig. |2^ we show a double dot geometry which 
is a candidate for the realization of the SU(3) x SU(2) 
model. In the leads the electrons can be described as free 
particles: 



loads 
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t=/lJ2 tcr 



€ CJ 



(5.36) 



where a = 11,12 refers to the two leads and a — ± is the 
electron spin. _ 
The Hamiltonian of the dots can be written asE£l 



ffdots -EE' 



a=dl,d 2 err 



{Qi-V^f , {Q 2 -V 2 C 2 f , Q X Q 2 



2C E ,i 



2C- 



S.2 



, (5.37) 



where a = dl, d2 is the index of the two dots, Qi and Q 2 
are the charges of them, Vi and V 2 denote the applied 
gate voltages in the Figure, rand the C's denote different 
capacitances of the system.E3 With a suitable choice of 
the gate voltages one can achieve that the ground state 
of the dots becomes three fold degenerate corresponding 
to the states (3 := (Qi,Q 2 ) = (0,0), (0,e) and (e,0). 
Then the tunneling processes among the leads and the 
dots result in the simultaneous flips in the electrons or- 
bital quantum number a — {11, 12, dl, d2} and the charge 
variables which now take over the role of the orbital 
index of the NLS. Since the tunneling is diagonal in the 
electron's spin we have an additional SU (2) degeneracy 
in the spin of the electrons. In sum, this system is a good 
candidate for the realization of the SU(3) x SU (2) model, 
where the SU (3) fixed point symmetry corresponds to 
the three fold degenerate ground states of the dots. 



V, 




dotl 
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FIG. 26. A mesoscopic double dot system, candidate for 
the SU(3) x SU(2) model. 

Similarly to the case of the two-channel Kondo 
modeljEa another possible realization of the SU(N) x 



VI. CONCLUSIONS 

We have studied the Multichannel Coqblin-Schrieffer 
model (MCCS) and its relation to the N-level system 
model (NLS). The properties of the MCCS model de- 
pend on both the spin and flavor symmetries, SU(N) 
and SU(f), as well as on the spin of the impurity. 

We have performed both analytical and numerical 
studies of the model. As with the multichannel Kondo 
model, there are three different classes of fixed points 
depending on the spin of the impurity, /x. The under- 
screened and completely screened fixed points (fj, > f 
and /i = /, respectively) have qualitative similar behav- 
ior to the analogous multichannel counterparts (N = 2). 

There are overscreened fixed points. They display 
Non-Fermi liquid behavior: they have associated anoma- 
lous residual entropy and anomalous exponents in the 
low-temperature expansion of quantities like the specific 
heat and the magnetic susceptibility. 

For an impurity with spin in the fundamental repre- 
sentation of SU(N), the residual entropy, S l , is only a 
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function of N+f and | log(iV//)|. Hence, there are differ- 
ent fixed points with the same value of S l . The exception 
corresponds to N = f, which yields the largest value of 
the residual entropy for fixed N + f. 

The low-temperature thermodynamics are determined 
by the value of the ratio 7 = f/N alone, for any /i < /. 
When / 7^ N we have 



— If-, x c ^ n - + 



T 



T 



(6.1) 



which diverge for 7 > 1, but remain finite for 7 < 1. 
When N = f, the power is replaced by a logarithm, as 
i n th e two-channel Kondo model. The constant terms in 
( |6.l[ ) are always present and are the dominant contribu- 
tions in the completely screened case, [i = /, and when 
/ < N. The Wilson ratio in such cases is given by 



R = 



N(N + f) 
N 2 - 1 



Channel anisotropy is a relevant perturbation. As the 
channel symmetry is reduced from SU(f) to SU(f'), the 
entropy is quenched since N+f decreases and | log(iV//)| 
increases. Likewise, a system with channel anisotropy 
might behave like a / > N system at intermediate tem- 
peratures and flow at low-T to a / < N system. 

Then we turned to the comparison of the MCCS model 
to the N-level system model (NLS), describing a heavy 
particle tunneling between N different positions and in- 
teracting with the conduction electrons. We have shown 
that the low-energy fixed point of the NLS model is 
just the SU(N) x SU(f) MCCS model. Performing a 
1/ f study of the NLS model we have analyzed the op- 
erator content of this low-energy fixed point, and the 
scaling properties of different physical quantities in the 
N < f limit. We have shown in this limit that while 
the operator content of the NLS model is different from 
that of the MCCS model, apart from some subtle differ- 
ences, the low-energy properties of the two models are 
the same. Especially, comparison with the exact results 
obtained in±he first part of the paper and with the NCA 
calculations^ show that the susceptibility, the residual 
entropy and the resistivity of the two models behave in 
the same way, and for reasonable physical parameters 
even the scaling of the specific heat is properly described 
by the MCCS model. 

Finally, we discussed some possible physical realiza- 
tions of the SU(N) x SU(f) models. First we discussed 
the case of tunneling interstitials in multicomponent met- 
als such as Pb\- x Ge x Te compounds. We pointed out 
that the low concentration of interstitials is essential to 
avoid strong inter-impurity interactions and keep the di- 
agonal elements of the self-energy A ab small. Secondly, 
we suggested a double quantum dot structure that could 
give an ideal realization for the SU(3) x SU(2) model. 
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